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Abstract. In a recent paper (Straus et al. 2008) we determined the energy 
flux of internal gravity waves in the lower solar atmosphere using a combination 
of 3D numerical simulations and observations obtained with the IBIS instrument 
operated at the Dunn Solar Telescope and the Michelson Doppler Imager (MDI) 
on SOHO. In this paper we extend these studies using coordinated observa- 
tions from SOT/NFI and SOT/SP on Hinode and MDI. The new measurements 
confirm that gravity waves are the dominant phenomenon in the quiet mid- 
dle/upper photosphere and that they transport more mechanical energy than 
the high-frequency (> SmHz) acoustic waves, even though we find an acousti c 
fiux 3-5 times larger than the upper limit estimate of Possum & Carlsson ( 2005 1 . 
It therefore appears justified to reconsider the significance of (non-MjHD waves 
for the energy balance of the solar chromosphere. 



Even sixty years after the two pioneering papers by [Schwarzschild ( 1948 ) 
and Biermann (1948), who independently proposed high frequency acoustic 
waves as a key agent, the origin and mechanism of the heating of stellar at- 
mospheres and coronae remains one of the great unsolved problems in solar 
and stellar physics. Recent studies by Fossum h Carlsson (2005, 2006), here- 



after F&C, and Carlsson et al. (2007) conclude that high frequency waves are 



not sufficient to heat the solar chromosphere. Others (e.g. Cuntz et al. 2007 
Wedemeyer-Bohm et al. 20071 Kalkofen 2007) question these results and argue 



for high-frequency waves to play an important role. 

Stably stratified atmospheres can support and propagate not only acoustic 
waves but also internal gravity waves. These have been proposed as an agent for 
the mechanical heating of stellar atmospheres more than 25 years ago (Mihalas &; 
Toomre |l981 ). However, the difficulties associated with directly observing them 
in our closest star, the Sun, have resulted in their neglect in this matter. Re- 
cently, by combining high quality 2D observations and 3D numerical simulations, 
Straus et al. (2008), hereafter SFJ, were able to determine the height dependence 



of the energy flux of internal gravity waves in the lower solar atmosphere and 
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found that at the base of the Sun's chromosphere it is around 5kWm~^, which 
is comparable to the radiative losses of the entire chromosphere. 

Here we revisit the role of both internal gravity waves as well as high fre- 
quency acoustic waves in the dynamics and energetics of the Sun's atmosphere 
using high quality Dopplcr velocity measurements obtained with SOT/NFI and 
SOT/SP on Hinode and with MDl on SOHO. 

The study of wave propagation characteristics requires simultaneous mea- 
surements in at least two different heights in the solar atmosphere. Gravity waves 
have long periods and small spatial scales, and a proper horizontal wavenumber 
separation is essential for their study. They therefore require both high spatial 
resolution and a large field-of-view, which limits the achievable cycle time of 
current instrumentation. For high frequency waves, on the other hand, a short 
cycle time is the key requirement. We therefore developed different observing 
programs, one targeted at gravity waves, the other at high frequency waves. 

To study gravity waves, we combined SOT/NFI filtergrams in Mg b2 with 
simultaneous MDI high-resolution observations in the Ni line at 6768 A. The 
data were obtained on 2007 October 20. Both time series extend over 12 hours 
with a cadence of 40 seconds and 60 seconds and a spatial resolution of 0'.'32 / pixel 
and O'.'61/pixel for SOT/NFI and MDI, respectively. The NFI filtergrams were 
interpolated to the times of the MDI observations using cubic splines. Dopp- 
lergrams were determined from differencing the blue and red wing filtergrams. 
They were calibrated by comparing the cumulative probability distributions of 
the Dopplergrams made with the normal filter settings, with those made with 
settings shifted by 20 mA to the blue. The resulting Dopplergrams were coregis- 
tered and mapped onto the lower resolution MDI data. Both data cubes v{x, y, t) 
were subjected to a 3D Fourier transform yielding complex Fourier coefficients 
v{kx,ky,v). Fig. 1 shows the phase difference spectrum ^(/c/i, z^) between the 
Mg b2 and Ni 6768 Doppler signals after azimuthal integration of the complex 
crosspower VMg{kx, ky, v) v^-^{kx, ky, u) over k"^ = k'^ + ky. It beautifully confirms 
the presence of internal gravity waves in the upper photosphere (transition to 
green-yellow-red in the gravity wave regime), and, thanks to the high wavenum- 
ber and frequency resolution, reveals interesting ridge/interridge structures in 
the region of p-modcs, which will be discussed in another paper. For a height 
difference of 600 km between the two signal forming layers (derived from the 
slope of the linear phase delay in the high-frequency part of the phase spec- 
trum), the phase spectrum yields the vertical phase velocity fph,z(fch; '^)) which 
in turn yields the vertical group velocity Vg^z = —Vph,z sin^ (9) {6 is the angle 
between the oscillation in the wave and the vertical, for details see SFJ). With 
the Mg b2 power spectrum and the VAL density of 5.97 • 10~^ kg m~^ at 720 km, 
the vertical energy flux carried by gravity waves at that height can be estimated 
as po{v'^)vgz{kh, ~ 2kWm~^, which is consistent with the results of SFJ (Fig. 
2). 

For the study of high frequency waves, we used two data sets obtained 
simultaneously with the spectropolarimeter (SP) and narrowband filter (NFI). 
For this short paper, we restrict ourselves to the discussion of the results from 
a data set obtained on 2007/08/13, and here on the results from the Fei 6301 A 
line, the stronger one of the two Fe lines recorded with SP. The data set extends 
over 237 minutes at a cadence of 16 seconds. There are 1024 pixels along the 
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Figure 1. V-V phase dif- 
ference spectrum ^{kh^ i^) be- 
tween Mg b2 and Ni i 6768 A. 
Note the area of downward 
phase propagation (green- 
yellow-red) in the region 
of propagating atmospheric 
gravity waves in the k — cu- 
diagram. 
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Figure 2. Flux of atmospheric gravity waves (red) and h igh-frequency 



acoust ic waves (green) in the solar atmosphere, as found by Straus et al. 
(2008) in numerical simulations (lines) and observations (symbols), 'i'he tlux 
has been integrated over the regions in the kh — v diagram shown in the insert. 
The new results obtained with Hinode in this work are added in bold. 
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0.01 Figure 3. Flux of high- 

frequency acoustic waves from 
observations with Hinode in 
the Fe 6301 line, before 
(black) and after (blue) cor- 
rection for the line formation 
MTF (red). The flux given 
in F&C is shown in green for 
comparison; see text. 
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slit with a pixel size of (/.'IG. The Doppler signal of this line is formed at 
around 400 km, i.e. close to the UV continuum analyzed by F&C. A critically 
important factor in determining the energy flux of high frequency waves is the 
line formation modulation transfer function (MTF; Keil &; Marmolino, 1986). 
To better estimate this effect, we model a spectrum of monochromatic waves in 
a stratified atmosphere with frequencies in the range 6 to 100 mHz and constant 
amplitude. With the synthesis of the Fe 6301 line profile the MTF can be 
determined as the ratio of the velocity power at the formation height of the 
line and the power of the Dopplershift velocity determined in the line core. 
The results are shown in Fig. 3, with the MTF in red, the measurements of 
the acoustic flux from F&C in green, the uncorrected flux determined from 
the power spectrum in black, and three MTF-corrected flux curves (assuming 
different S/N levels in the power tail) in blue. Note the steep fall-off of the 
MTF of this typical photospheric line by two orders of magnitude in the range 
5 < < 30 mHz. Depending on the noise level, taking into account MTF effects 
could even lead to a turnaround of power spectra at high frequencies. Using the 
thicker, middle blue curve (assumed S/N = 0.17), we estimate the flux of waves 
with frequencies between 5 and 31 mHz at 400 km height to be 1600 Wm~^ 
(550Wm~2 for 10 < z/ < 31 mHz, and 13Wm~2mHz~^ for v > 31 mHz), which 
is significantly larger than the estimate of F&C (< 438 Wm^^ for waves between 
5 and 50 mHz). 

We feel that it is justified to reconsider the significance of the energy flux of 
(non-M)HD waves, for four reasons: (1) In the past, the contribution of internal- 
gravity waves has been underappreciated. They turn out to carry enough energy 
to balance the radiative losses of the chromosphere. (2) With Hinode data we 
find an acoustic flux 3-5 times larger than F&C. We believe that their estimate 
of the upper limit is underestimated. (3) The effects of the line formation MTF 
seem to have been underestimated in the past. (4) The propagation behaviour 
of high frequency acoustic waves in the "real" solar atmosphere appears to be 
more complicated than commonly assumed. 
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